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Abstract 

A state coupling between the hydrostatic (volumetric) and deviatoric 
parts of the free energy is introduced in a damage mechanics model rele- 
vant for the quasi-brittle materials. It is shown that it describes the large 
dilatancy of concrete under compression and the different localization an- 
gles and damage levels in tension and compression. A simple isotropic 
description is used, although similar ideas can be extended to anisotropic 
damage. The model is identified with respect to tensile and compression 
tests and validated on bi-compression and bi-tension. Fully written in 
three dimensions under the framework of thermodynamics of irreversible 
processes, it allows further developments within a finite element code. 

1 Introduction 

At the ultimate state of damage, quasi-brittle materials can be regarded as gran- 
ular materials whose behavior, due to grain interlock, is strongly dependent on 
the confinement. Classical yield criteria (Mohr-Coulomb, Drucker-Prager,...) 
use the confinement (the hydrostatic strain) as a reinforcement factor. But 
there is experimental evidence that elasticity of powder materials, thus their 
free energy, also depend on it [1] . That point is generally not taken into account 
in damage models. As a consequence, the large dilantancy observed at these 
states, for example on concretes [2, 3], is poorly described by classical damage 
models as apparent Poisson's ratio cannot exceed 0.5, the limit for linear elas- 
ticity, while experiments exhibit greater values. Another specificity of concrete 
behavior is the very different crack angles observed in tension and compression. 
In most damage models, the Hadamart and Rice criterion [4, 5] leads to the same 
localization angle in both cases. Furthermore, for a concrete specimen, even af- 
ter a rupture in tension, some carrying capacity in compression remains: this 
implies a damage level at the onset of localization much lower in tension than 
in compression. The present model describes these effects that are generally 
missed by most of damage models. 

This work, in continuity with [6] , constitutes an attempt to use a non convex 
potential in the field of damage mechanics. For this reason, a simple isotropic 
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damage law is considered. As in the Kelvin's approach of elasticity [7], plasticity 
theory, soil mechanics and some damage models [8] , the isotropic and deviatoric 
decomposition is used. The retained yield criterion [9] is smooth and convex. 
The model is identified with respect to the well known uniaxial and multiaxial 
testings of [3]. 

2 Constitutive law 

Helmholtz free energy. The damage level is described by the scalar variable 
d that ranges from for sound material to 1 for fully damaged material [10]. 

The present model is an isotropic one: the hydrostatic and deviatoric partitions 
of the stress G = g^+gf and the strain e = e*^ + e*^ are used (details are given in 
Appendix A). The state variables (e'*,£'^,d) describe the material's state and, 
with respect to the generalized standard material framework [11], the associated 
thermodynamic forces are respectively (a'^, a*^, Y) where Y is the energy release 
rate density. The proposed free energy ^ is: 

2p^ie'\i', d) = 3Ke}' : e*^ + 2^ (l - d{l + 2ipe'^)) e'^ : ef, (1) 

where e'' = ti{e)/V^ (Appendix A). The new constant introduced is if, the 
other ones are the mass density p and the bulk and shear moduli, respectively 
K and /i. The role of cp will be detailed further but it can be already seen that 
it introduces a cubic term in the free energy and that setting ip = leads to a 
simple damage model in which only the deviatoric part is aff'ected by damage. 
This choice has been made in order to get rid of the unilateral effect of damage 
on bulk modulus (that comes from crack opening and closure) because the use 
of positive parts of the strain tensor induces difficulties associated with non 
regular free energy [12]. In Eq. 1 tensors e}^ and e'^ can be replaced by their 
algebraic values and e'' (Appendix A). 

Stress to strain relationship. The hydrostatic and deviatoric stresses, as 
thermodynamic forces, are obtained by differentiation of the free energy with 
respect to e}^ and e'^: 

= 3Ks^ - ^ij,ipd{e'^fL, (2) 
v3 

a'^ = 2iJ,{l-d{l + 2ipe^))s'^. (3) 

Setting d = leads to recover the linear isotropic elasticity law in the Kelvin's 
decomposition form [7]. The deviatoric stress and strain remain coUinear to- 
gether (and coUinear to the unitary tensor D^, see Appendix A) . Once projected 
on the orthogonal tensor base (X,:D), the previous expression becomes: 

= 3Ks^-2fi(pd{e'^f, (4) 
a"^ = 2i^{l-d-2ipde^)e'^. (5) 

Let us suppose an imposed deviatoric strain while the hydrostatic strain remains 
equal to zero (isochoric transformation) then a confining pressure (cr^ < 0) is 
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necessary to keep the volume unchanged: 



C7*>(£*' = 0, s") = -2fi<pd{ey. (6) 

Let us suppose now that the deviatoric strain is joined with an hydrostatic stress 
imposed equal to zero then an induced dilatation e*^ > arises: 

,v=o,.-)=>f&;£^ (T) 

These two aspects are relevant to the dilatancy effect in concrete-like materials. 
The micro-mechanical point of view associated to these effects is the surmount- 
ing of concrete particles, inducing voids creation, that arises under irreversible 
shearing [13, 14]. 

Elastic tangent modulus. The tangent modulus is defined as H = dq_/de. 
In case of elastic transformation, d remains fixed, then H reduces to M° = da/de. 
Prom the stress to strain relations (2, 3), we have: 

HO = ^^~^^ 7g) J + 2/2I-a(Jg)£^+£'^g) J), (8) 
2/i = ^ = 2^1 - d) - 4/.(pde^ (9) 

In this expression I is the fourth order identity whose expression in index form 
is [Ijijki = (^ikf^ji + 5ii5jk)/2, the symbol (8) refers to the tensor product and 
(from Eq. 5) jl represents the apparent shear modulus. Then If has the index 
symmetries of an elasticity tensor: [H*']jjki = [IHI°]kiij = [H^Jijik- 

Inverse stress to strain relationship. When d ^ 0, (the case d = is 
straightforward), replacing in Eq. 5 by its value (Eq. 4) gives, first, s'^ as the 
solutions of a third order equation and, second, the relation between and £**: 

{s'^f+ps'^ + q = 0, 
2ipda^ - SK{1 - d) 3Ka'^ 



, _ a^ + 2^^^d{e<'y 

3K ■ ^ ' 

The discriminant A = (p/3)^ + (5/2)^ leads to analyse from one to three real 
solutions. This non univocal inverse relationship would not be used in finite 
element calculation, but we shall have to consider these three possibilities for 
the semi-analytical resolutions in Sec. 5 and 6. 



3 Free energy analysis 

The proposed Helmholtz non convex free energy exhibits complex behavior that 
has to be carefully studied. Fig. 1 shows the contour plot of the free energy 
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defined by Eq. 1. The transformation APCZ is tlie compression studied in 
Sec. 5 and the damage value retained for the drawing is d{C), the maximum 
reached during this transformation. 




Figure 1: Iso- values of the free energy in the Drucker-Prager plane for d=0.64 



Positivity of the free energy. It signifies that the material cannot restore 
more energy than has been stored inside. In linear elasticity, it is associated with 
the positivcness of the Kelvin's moduli (3X, 2/i) and infers classic bounds for 
the Poisson's ratio (— 1 < ^ 0.5). In the general case, from Eq. 1 we have: 



<0 4^e'^ > ^-^e^. (12) 

This condition corresponds to the forbidden domain inside the curve p'^ = 
in Fig. 1. It only exists when p. < 0. It will be shown in Sec. 3 that it cannot 
be reached during a transformation because localization occurs at least when 

fl = 0. 

Particular lines. Lines of interest are the loci of a'^ = and = 0. From 
Eq. 5 we have: 

= ^ = 0. (13) 

Prom Eq. 9, this corresponds to p, = 0. On the right side of this line on Fig. 1, 
p < 0. The second solution for ct'* = is more classically = 0. The locus of 
= (Fig. 1) is given by Eq. 4: 



ed^J|^^ah = 0. (14) 
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Between the deviatoric axis = and this hne, > and cr^ < 0. This 
represents the dilatancy effect of the modeh the deviatoric strain induces dila- 
tancy even if the material is under (moderated) pressure. The saddle point H 
is defined as the intersection of the two lines cr^ = and cr*^ = 0, we have: 

As g_{H) — 0, it corresponds to an unstable free stress state with £{H) ^ 0. 
Finally, it can be easily shown that the point N, corresponding to the minimum 
deviatoric strain for = is such as: £^(N) = 26*^ (H), e^{n) = V2e'^(H). 

Positivity of the energy release rate. The energy release rate Y is 
associated with the damage d: 

Y = -^=n{l + 2^e'^){e''f. (16) 

The thermal dissipation is qf = Yd. As damage d cannot physically decrease, 
its positiveness implies that if Y is negative, the evolution of the damage must 
stop: 

< — ^d = 0. (17) 

In the example of Sec. 5, = 160 then this condition is e*^ < —3.12 10~^. That 
region, denoted as d = on Fig. 1 corresponds to high confinements out of 
most practical applications (far to be reached in the presented examples); it 
is strongly due to the initial choice of damage acting only on the hydrostatic 
part. One can remark that some experiments tend to confirm the existence of 
a damage locking at high confinements [15]. 

Localization. The Hadamart and Rice criterion [4, 5] of localization au- 
thorizes the existence of a localization plane (a concentration of large strains) 
orthogonal to the unit vector n as soon as: 

det(n.H.n) < 0. (18) 

Although we still consider here elastic transformations, we assume this local- 
ization to correspond to the apparition of a macroscopic crack. Once appeared, 
the damage model ceases to apply as the body is split into pieces: it infers a 
restriction to the domain of definition of the model. Among possible vectors n, 
we consider nj, an eigenvector of the strain deviator tensor, then ^.n\ = efni, 
where is the corresponding eigenvalue of e'^ (a principal deviatoric strain). 
As e = e° in this case, Eq. 8 give: 

niM°.ni = ( ^^^^'^ - 2Q;ei^ ni ® ni + /iJ. (19) 



The determinant of this expression is: 

det(ni.H°.ni) = fi^ ( ^E±^ _ 2aef ] . (20) 
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This shows that /i = is a sufRcicnt condition of locahzation. The domain 
(grayed) on the right side of this hne in Fig. 1 will not be concerned by the 
possible evolutions. Prom Eq. 12, it always contains the domain ^ 0. From 
Eq. 63, ef ^ ^ 0, considering that K, jl and a are positive and using 

Eq. 9 and 10, we obtain a second sufficient condition for localization: 



2,K + 4/x(l - d) 



(21) 



The corresponding line is denoted as ry on the Fig. 1. The localization occurs 
at least when a transformation reaches the lines or /i = 0; the domain above 
them (grayed) is not concerned by possible evolutions. 

Other remarks about stability. The current line w passes by the saddle 
point H. From its definition e^ja^ = e^ju^ comes: 



2^(1 - d)e'^eh - 4/i(^de'^e'^eh = 'iKe^e'^ - 2tiipd{e'^) 



(22) 



This expression has no simple analytic solution and ui has been drawn on Fig. 1 
with a steepest slope algorithm. Below u, the stress tends to bring back the 
system to the stable original state 0; above u it tends to make the system reach 
the domain above the line rj where the localization occurs. Then the domain 
between oj and r] can be considered as unstable, leading to an instantaneous 
evolution towards the lines t] and jl = where the localization occurs. 

4 Dissipative behaviour 

Yield surface. This yield surface is based on the elastic criterion proposed 
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Figure 2: Initial yield surface (MPa) 



by [9]. Its shape, smooth and convex, can be regarded as a softened approxi- 
mation of the Von Mises and Rankine criterions. It is in good agreement with 
the elastic limit identified by [3] in biaxial testings (Fig. 2). The first member 
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corresponds to the Von Miscs expression; the second one, due to a tensor expo- 
nential, growths quickly for positive values of any principal stress. The material 
constants are try, that principally rules the limit stress in compression and tro, 
that strongly influences the tension to compression stress ratio. The fimctions 
g{d) and h{d), with g{0) = h{0) = 1, are hardening functions (whose detailed 
expressions are given by Eq. 47 and 48). 



f{QL,d) = 



+ 



9{d) h{d) 



- V3 I -cr- 



(23) 



Naming (cji, an, am) the eigenvalues of the stress tensor and considering sound 
material (d = 0), this criterion writes: 



/((Tl,(7ii,£7iii,rf = 0) = 



(cjii - tTiii)2 + (ctiii - cri)2 + (cti - any 



'2ai\ , /2aii 
CTo I xjexp [ — 1 +exp I — — 



-exp(^)-y3l 



The identification of ((Jq, Cy) is done with respect to the measured elastic limit 
stresses a'^ = —9.98 MPa in pure compression and a* = 1.29 MPa in pure 
tension. The difference f{a\0, 0, 0) - /(a^ 0, 0, 0) gives: 



aoWexp 



2crt 



(JoWexp 



(25) 



This equation accepts a numeric resolution given by the intersection of the two 
functions of ao that correspond to left and right members. This intersection 
is unique because each one is monotonic. The value of ay is obtained from 
/(cr*, 0, 0, 0) = or /(cr'^, 0, 0, 0) = 0. Numerical values of (ao, Cy) are reported 
in table (1). 

Tangent modulus and evolution equation. In case of dissipative trans- 
formation, the consistency equation df{a, d) = and the differential of a{s, d) 
give the damage evolution dd: 



dd=- 



dd da 



dd 



da 



¥0 . 



: de. 



and the expression of the tangent modulus H = da/de: 



1 + 



dfY^ 9a df 



dd 



dd da 



(26) 



(27) 



The inverse of the fourth rank tensor can easily be obtained in a convenient 
tensorial base [7] as the inverse of a 6x6 square matrix. The involved derivatives 
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are obtained without difficulty from Eq. 2, 3 and 23: 



dl 
dd 

dl 
da 

dd 



\a^\\g'{d) aoh'id) 



+ 



K^id) 
1 exp (2^/(70 ) 



g_ 

exp ( — 

CTO 



9{d)\\a^\\ h{d) ||exp(a/ao)||' 
-2M£'')'^-2M1 + 2^£V- 



V3 



(28) 
(29) 
(30) 



The elasticity equations (2, 3), the yield condition (23) and the evolution equa- 
tions (26 to 30) constitute the set that is required for the use of the model within 
a FEM code. 



5 Tension and compression simulation 



stress to strain relationship. The stress tensor expresses as ct = ae\ e\. 
Due to the isotropy of both the material and the model, the strain tensor writes: 



£00 

-v£ 
-ve 



(31) 



where e is the uni-dimensional strain and v the apparent Poisson's ratio. Using 
the projection equations (57, 58) and denoting D = £'^/||£'*|| (see Appendix A), 
we have: 



= sign(£(l + v)) 



£^ = 



1 -2z> 



^/3 



-£, £ 



2/\/6 
-1/V6 

-i/Ve 

^^|£(1 + Z^)|. 



(32) 



(33) 



Using these equations in Eq. 4 and 5 leads to two expressions for the stress to 
strain relationship: 



a = ZK{l-2i>)E~^{l + ufE\ 
v3 

<T = 2n(^-d-2ipd^^^^£^{l + v)e. 
Eliminating a between them gives a second order equation in v: 



(34) 
(35) 



2 12 uipd 
^ 7=- ' 



-e + v 



2ii{l-d) + 6K + 



12 jjLipde 



+ 2n{l ~d)-3K = 0. (36) 



Eq. 34 or 35 gives from zero to two solutions a for a given strain s (these 
solutions correspond to different values of i/, then to different strain tensors e). 
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Dissipative trEinsformations equations. The yield function (23) be- 
comes: 



CTO 



^9{d) h{d) 



I exp 



+2-^/3j -ay. (37) 



As in the general 3D case, from the differential of a{e,d) and the consistency 
equation d/* (a, d) = is obtained the stress to strain differential relationship: 



da 



df* df*da\ df*da 



dd 



+ 



Prom Eq. 37 we have: 



da 
dd 



2 sign(o-) 

3 gid) 



da dd 



exp 



dd ds 



ds. 



+ 



h{d)^2 + e^p[^) 

2\a\g'{d) aoh'{d) ' ' 



3 g-'id) h^d) 



2 + exp 



V3' 



(38) 



(39) 



(40) 



From the stress to strain relationship (Eq. 34,35 and 36), the following deriva- 
tives, in which da/de is the tangent Young modulus and A the discriminant (in 
Eq. 36), are obtained: 



da_ 

de 
da 



2/is{l + v) dd 



ZK{1 - 2u) - -^IJ-fdil + i>fe ■ 

^|(2.|-.l)a..,.- 



'iK+^,jL^d{l + v)e]2e^, (41) 



1 + 2^?- 



^/3 

1 - 2z> 



de' 



di 



2(1 + + 2£;^4)2) 



-{l + 2u)± 



-2/i(l -d) + 12iiifd£/\/Z 

• 



(43) 



Localization in tension and compression. In the tension or compression 
state, the axis 1 is the symmetry axis. As axis 2 and 3 have an equivalent role, 
we can define with no restriction n in the plane [61,6*2] as: n = cos(^)ei + 
sin(0)e2. For dissipative transformations, the general equations of Sec. 4 allow 
to compute det(n.H.n) numerically by searching its minimum value with respect 
to 6. The non convex potential authorizes possible localization during elastic 
transformations. Prom Eq. 8 we have: 



n.H°.n = 



h+(a + 4c) cos2 (61) (a + c) sin(6') cos(6l) 
(a + c) sin(6l) cos\e) h+{a- 2c) sin^ (6*) 
6 



(44) 



a = i^+|(l-d)-^M(^rf(l-2i>)£, 



b = jiil—d) —fnpd{l—2i/)e, 

v3 



c = — 



3\/3 
(45) 



IJ,(pd{l+i')e. 
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Identification. The identification of tlie Young Modulus E = 31.2 GPa and 
the Poisson's ratio v = 0.162 being very classic, they will not be detailed; the 
corresponding bulk and shear moduli {K,n) are given in table (1). The proce- 
dure of identification of (cto, o-y) is given in Sec. 4. Eliminating d between Eq. 34 
and 35 gives the following expression that can be used with the experimental 
data in order to determine ^p: 

_ ^/3 a/e - '6K{1 - 2v) 

^ 2 2,Ke{l-2vy-2^ie{l + vY+a^>' ^ ' 

The hardening functions g{d) and h{d) are identified with respect to the trans- 
formation AC (see Fig. 3). Their structure has been chosen in order to have a 
smooth transition at the yielding point A. 

g{d) = l+5ov^, (47) 
h{d) = l-d + ho{y/d-d). (48) 



Table 1: Material constants 



K n 




50 ho 


15.4 13.4 
(GPa) (GPa) 


8.004 0.4551 160 
(MPa) (MPa) 


3.00 160 



Uniaxial compression simulation. The set of equations (31 to 43) defines 

the strain driven compression and tension curves (Fig. 3 and 4), allowing the 
comparison with Kupfer's data (circles). The compression starts by the linear 
elastic transformation OA, with = 0. The material begins to yield at point A 
where (t(A)= a'^ . The transformation ABC corresponds to the damage growth, 
up to point C where d{C)= 0.64. The localization criterion is never reached 
during ABC. At point P, z> = 0.5: the material begins to dilate. Unloadings 
such as OBO for d = 0.5 and OCO for d = 0.64 reveal weak nonlinear elasticity: 
it can be noticed that the curvature of OC (the tangent modulus decreases as the 
load increases) is qualitatively in good agreement with the experiments of [16] for 
example. Of course hysteresis loops or permanent strains are not described by 
the present simple form of the model that does not takes into account plasticity 
effects. Point C corresponds to a horizontal slope for the transformation. Then 
da = and, as the value of df /dd 7^ at this point, Eq. 38 shows that da/de = 
0: the elastic transformation OCZ has a also an horizontal tangent at point C. 
After C, da/de < 0, the transformation CZ is a non linear elastic softening. The 
damage remains fixed at d{C). From some point just before Z, an unloading 
follows the unrealistic elastic transformation ZCO (this aspect of the model 
remains to be enhanced). At point Z the localization condition relative to 
the elastic transformations (Eq. 45), is reached and this point is close to the 
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£11-2 2 £22 = £33 (%o) 

Figure 3: Uniaxial compression simulation 



experimental breakdown point identified by Kupfer (although this value has 
not been used during the identification process). The continuation of the elastic 
transformation after Z is, for information, drawn on Fig. 1; as and e*^ are 
dependent upon s by Eq. 33, they become simultaneously equal to zero at the 
saddle point H. The localization arises at point Z for the angle 9'^ = 56 degrees: 
the crack plane is rather close to the load axis ei, as classically observed in 
experiments. The apparent Poisson's ratio v reaches 0.70 at the point Z and 
this large value, that cannot be reached by a classic damage model, is responsible 
of the good description of the transverse strain. 

Uniaxial tension simulation. The tension curve (Fig. 4) is also referred 
to points A,B,C',Z whose roles are similar to the ones used in compression. 
The yielding occurs at point A, at the identified value cr* that is, according to 
Kupfer and confirmed by [17], less than the half of the peak stress: this result 
is different from the classic vision of a fragile behavior in tension. At point C, 
df /dd = then (from Eq. 38) the slope is horizontal. It is not reached because 
the localization occurs before, at point Z, during the dissipative transformation. 
If [17] describes a post-peak evolution, Kupfer does not, neither this model with 
the retained constants. Again, although not used in the identification process, 
the stress (j(Z) is close to the experimental value. The angle of localization at 
point Z is 0* = 27 degrees, leading to a crack plane rather close to be orthogonal 
to the tensile axis e\, as commonly observed (moreover most of damage models 
suppose a crack at degree). The lateral strain, due to low damage values, 
remains very small, consistently with experiments. The damage value at the 
localization point (i(Z)= 0.22 is weak: although isotropic, this model allows to 
have a rupture in tension and to keep a loading capacity in compression. 
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d(Z) = 0.22 

= 27= 



T 



£22 = 633 0.05 0.1 eii %o 

Figure 4: Uniaxial tension simulation 



6 Equi-biaxial simulations 

In this section we consider an equi-biaxial loading (strain driven) in order to 

validate the model with the constants identified in tension and compression. 
The stress tensor is a = a{e2 (81 62 + 63 (81 63). Due to the initial isotropy of both 
the material and the model, the strain tensor writes: 



-ue 

e 
£ 



(49) 



with £ < 0. The projection equations (57, 58) give: 
sign(£(l + />)) 



D 



-2/V6 
1/V6 







1/V6 J 



%/3 ^' 



\e{l + i>)\. 



(50) 



(51) 



Using these equations in Eq. 4 and 5 gives two expressions for the stress to 
strain relationship: 



- = f (2-%-^(l + .)V, 
a = 2/i ^1 - - 2<fid^^^e^ (1 + i>)e. 



(52) 
(53) 
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Eliminating a gives the following second order equation which accepts from zero 
to two solutions in v: 

i''^2Vllj.^de + V ^2/i(l -d) + + 2/x(l -d)-ZK- 2\/^n^de = 0. (54) 

Other calculi are similar to the previous case. 

Equi-biaxial compression simulation. Directly dependent upon the 




£22 -£33 -2 2 4 £11 (%o) 

Figure 5: Equi-biaxial compression simulation 



yield surface, the yield point A (Fig. 5) is well depicted. The damage increases 
during transformation AZ up to d{Z)= 0.70, this value being greater than the 
maximal one obtained in uniaxial compression. At point Z, (which is before the 
point C where df /dd = 0) the localization condition for dissipative transforma- 
tions (Eq. 18) is reached. The localization angle is 9'^'^ = 40 degrees, consistently 
with Kupfer's post mortem picture. The conical crack shape in Fig. 5 has been 
chosen in order to respect the symmetry axis of the problem. Without contra- 
diction, the observed structure is pyramidal, this shape allowing the kinematics 
of the concrete pieces. The localization appears too quickly and the large dila- 
tancy is a little overestimated. Moreover, this post peak region differs according 
to different experimental methods and concretes: for example the high strength 
concrete tested by [18] exhibits no post peak. 

Equi-biaxial tension simulation. The part before the peak stress is in 
good agreement with the experiments for both axial and lateral strains (Fig. 6). 
Contrary to experiments, a dissipative post-peak evolution C'Z is described (it 
represents a minor problem compared to the elastic transformation CZ in the 
uniaxial compression case). The damage d{Z)= 0.67 at ultimate point reaches 
about the same value than in biaxial tension. The localization angle = 12.6 
degrees obtained in this case; involves cracks approximatively orthogonal to the 
load directions, very consistent with Kupfer's post mortem pictures. Again, the 
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-0.1 0.1 %o 

Figure 6: Equi-biaxial tension simulation 



conical shape is only the representation of the symmetry of the problem as the 
model only defines the angle of the crack, not its shape. 
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7 Conclusions 



This study shows that the simple coupling introduced between the hydrostatic 
and deviatoric parts of the strain in the free energy leads to a mimber of conse- 
quences, even if the primary effect, the non linear elasticity, slightly influences 
the shape of the elastic discharges. The important dilatancy of the concrete 
at ultimate stages is well depicted. The localization criterion is used as an in- 
dicator of the creation of a macroscopic crack, it applies at stress and strain 
levels that are quite consistent with experiments. The damage value predicted 
at localization in pure tension remains small, letting to the material some load 
bearing capacity, for example for further compression. The crack angle given by 
the localization condition is different in tension and compression and consistent 
with observations. 

The description of these effects constitutes an improvement with respect to 
most of available damage models of comparable complexity (five constants are 
used). They are particularly relevant when, as in earthquake engineering, the 
response of a concrete structure under severe conditions is considered. The 
global structure of the model and its embedding thermodynamic framework, 
opens the way to numerical implementation within a finite element code. 

The introduction of a plasticity formalism [19, 20] would help to describe 
missing efi'ects such as hysteresis loops and permanent strains and to correct 
some irrelevant post-peak responses. The hydrostatic part of the free energy 
could be affected by damage in order to extent the field of application of the 
model to confined states and avoid the damage locking condition present in the 
actual form. The three-dimensionnal stress to strain relationship is univocal, 
however the inverse one is not: interesting theoretical analysis can be done on 
stress driven loadings as they may involve multi phased states (this refers to the 
work of [21] in elasticity and [22, 23] in plasticity). 

Appendices 

A Tensorial formalism 

For any symmetric second order tensor A the hydrostatic and deviatoric parts 
are denoted respectively by ^ and . The hydrostatic part is obtained by 
the projection onto the normed hydrostatic tensor // VS (where / is the identity 
tensor). 



The symbol ":" denotes the double contraction i.e. A : ^ = ^ij^ji = tr(^). 
In this model, the deviatoric parts and (of the Cauchy stress q_ and the 
infinitesimal strain e) will remain coUinear but may be a priori of opposite 



A = £'+A'^ 



(55) 




(56) 
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directions. We choose to define this direction with respect to the normed tensor 
given by: 




The norm used is the cuclidean natural one i.e. = Ai^A^y The algebraic 
values (£'^,£'*,C7'^,C7'^) are defined as: 

£^ = e:^, (T^=a:^, e"^ = e- D, a'^=a:D (58) 
V 3 V 3 

One can remark that e'^ is positive but a'^ can a priori either be positive or 
negative. These projections are known under the names p and q where p = 
—a^/\/3 and q = ^/3/2\a'^\ in soils mechanics. 



B Useful property of the deviator 

The sorted eigenvalues of e'^ are (ef ,efj,£fjj). Prom Eq. 58, we have: 

{4f + {sfif + iefuf = {e^f (59) 
4 + 4 + 4u = (60) 
4 > 4i > 4i (61) 

This leads to the following expressions: 

As soon as ej\ evolves between its bounds [—e'^/\/6;e'^/V6] (given respectively 
by £ji = £iii and £ii = £i ), these functions are monotonic and: 

(63) 
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